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BESCHREIBUNG 

Magne1resonaiiz-Ix)kaUsierungsverfahren 

Die Erfindung betrifft ein Magnelxesonanzverfahren zur Lokalisierung von interven- 
tionellen Vorrichtungen, insbesondere in vivo, bei dem die interventionelle Vorrichtung 
5 eine Markierung tragt, die bei Magnetresonanzaufhahmen die Messsignale beeinflusst 
oder eigene Messsignale erzeugt 

Die Anwendung von Magnetresonanzverfahren (MR-Verfahren) bei medizinischen 
Eingriffen gewinnt zunehmend an Bedeutung. Zum einen zeichnet sich die MR- 
10 Bildgebung durch hervorragenden Weichteilkontrast und dutch beliebige Orientierung 
der Bildebenen arts, zum anderen wild eine Gesundheitsgefahrdung der Patienten und 
des Bedienungspersonals durch ionisierende Sttahlung, wie sie bei Rontgenverfahren 
verwendet wird, vermieden. 

1 5 Bei der Visualisierung und Lokalisierung von interventionellen Vorrichtungen zur 

Einfuhrung in den Korper eines Patienten, insbesondere Kafhetern, stellt sich jedoch das 
Problem, dass diese meist nicht unmittelbar beobachtbar sind. Wahrend bei 
Bildgebungsverfehren, die auf Verwendung von Rontgenstrahlung basieren, selbst sehr 
Heine metallische Drahte einen ausreichenden Bildkontrast zur Kathetervisualisierung 

20 hervorrufen, bewirken diese in der Magnetresonanzbildgebung lediglich eine 

unzureichende Signalreduzierung, da derartig kleine Objekte nur ein sehr kleines 
Wasservolumen verdrangen. Aus diesem Grand muss die Sichtbarkeit der 
interventionellen Vorrichtungen auf andere Weise erhoht werden, wozu verschiedene 
Verfahren entwickelt wurden. 

25 

Die in der Literatur beschriebenen Lokalisierungsverfehren werden in zwei Kategorien 
eingeteilt. Bei den aktiven Verfehren wird die interventionelle Vorrichtung mit einer 
Empfangsspule ausgestattet, so dass Signale aus der Umgebung der Vorrichtung iiber 
einen zusatzlichen Kanal empfangen werden konnen. Passive Methoden hingegen 
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visualisieren die interventionelle Vorrichtung im MR-Bild durch den Kontrast 
gegenuber dem umgebenden Gewebe. 

Im Bereich der aktiven Verfahren haben sich bislang zwei Kameterlokalisierungs- 
5 verfahren etabliert. Zum einen kann in die Katheterspitze eine kleine Empfangsspule 
eingebaut werden, die fiber ein Koaxialkabel durch den Katheter hindurch mit einem 
Empfangskanal verbunden wird (C.L. Dumoulin et aL, Magn. Reson. Med. 29, 41 1-415 
(1993)). Der groBe Vorteil dieses Verfahrens ist die MSglichkeit, durch Anlegen von 
Feldgradienten die Koordinaten der Katheterspitze aus Projektionen in die entsprechen- 
10 den Raumrichtungen zu bestimmen. Das Verfahren ist daruber hinaus kompatibel mit 
alien schnellen Bildgebungsverfahren und somit echtzeitfahig. 

Alternativ zur Verwendung einer Empfangsspule kann auch eine elongierte Antenne in 
den Katheter eingefUhrt werden, die dann MR-Signale entlang des Katheters empfangt. 
1 5 Auf diese Weise konnen auch Instrumente mit geringem Durchmesser wie Fuhrungs- 
drahte und neurologische Katheter sichtbar gemacht werden. Ein besonderes Anwen- 
dungsgebiet ist die intravaskulare Bildgebung. 

Nachteilig macht sich bei beiden Verfahren bemerkbar, dass die durch den Katheter 
20 hindurch zum Empfangskanal verlaufende Leitung fur HF-Anregungspulse unbeab- 
sichtigt als Antenne wirken kann. So wurde gezeigt, dass sich ein Fuhrungsdraht nach 
30 Sekunden einer Gradienten-Echo-Sequenz auf 74° C erhitzen kann. Die Resonanz- 
bedingungen in diesem Fall sind vielfaltig und in der klinischen Praxis schwer zu 
kontrolheren. 

25 

Auf der anderen Seite stehen die passiven Techniken, bei denen die Sichtbarkeit der 
Katheter auf bestimmte Weise erhoht wird. Eine MogUchkeit ist die Verwendung von 
Kontrastmitteln, wobei entweder Katheter verwendet werden, deren Volumina mit 
einem entsprechenden Mittel (Gd-DTPA) gefullt sind oder deren Hulle 
30 kontrastverstarkendbeschichtetist. 
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Ein anderer Ansatz besteht darin, durch Storung des statischen Magnetfeldes B 0 
Suszeptibilitatsartefakte im MR-Bild zu erzeugen. Konventionelle Polyemylenkatheter 
konnen dazu mit paramagnetischen Ringen (Dy 2 0 3 ) prapariert werden. Eine Arbeits- 
gruppe am UniversitatsHinikum der RWTH-Aachen bat ein alternatives Verfahren 
5 entwickelt, bei dem eine lokale Feldinhomogenitat durch eine Drahtscblaufe im 
Katheter hervorgerufen wird, die an eine externe Stromquelle angeschlossen ist (A. 
Glowinski et al., Magn. Reson. Med. 38, 253-258 (1997)). Auf diese Weise ist das 
Bildartefakt fiber die Quelle wahrend der Intervention regelbar. 

10 Bei diesen drei passiven Visuabsierungstechniken ist positiv zu bewerten, dass es 

moglich is^ die gesamte Katheterlange sichtbar zu macben und das die Methoden mit 
alien Bildgebungsverfehren kompatibel sind Nacbteilig ist, dass alle Verfahren 
vergleichsweise zeitintensiv und die Koordinaten der Katheterposition nicht direkt 
zuganglicb sind. Folglich ist eine automatisierte Katheterverfolgung nicht moglich. 

15 

GemaB einem weiteren von M. Burl, Magn. Reson. Med. 36, 491-493 (1996) sowie S. 
WeiB, Proc. ISMRM, 544 (2001) beschriebenen Lokalisierungsverfahren wird ein 
Katheter, der auch als OptiMa-Katheter bezeichnet wird, an seiner Spitze mit einem 
elektronisch isoUerten Schwingkreis ausgestattet, der auf Larmorfrequenz abgestimmt 

20 ist. Wenn ein Bi-HF-hnpuls gesendet wird, wird der Schwingkreis angeregt und 
verursacht eine lokale Resonanzuberhohung des Bi-Feldes, die den FUpwinkel und 
somit das Signal lokal erhoht. Der Schwingkreis kami optisch fiber eine Photodiode, die 
. mit BQlfe eines durch den Katheter verlaufenden Iichtleiters beleuchtet wird, verstimmt 
und die Signalverslfirkuhg somit ah- und ausgeschaltet werden. Durch Subtraktion eines 

25 An-/Aussignals wird der Signalhintergrund unterdrfickt. Die erhaltenen Messsignale bei 
aktivierter bzw. desaktivierter Markierung werden auch als An- bzw. Ausprojektion 
bezeichnet 



Dieses Verfahren zeichnet sich dadurch aus, dass die Katheterkoordinaten direkt 
30 zuganghch sind und die Technik mit alien Bildgebungsverfehren kompatibel ist. Auch 
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die Patientensicherheit ist gewahrleistet, weil ein durch den Katheter verlaufender 
Lichtleiter im Gegensatz zu einem elektrischen Leiter nicht als Antenne wirken kann, 
die sich bei Einwirken von HF-Pulsen stark aufheizt SchlieBlich ist das Verfahren auch 
echtzeitfahig. 

5 

Nachteilig macht sich jedoch auch bei diesem Stand der Technik bemerkbar, dass die 
Detektion der interventionellen Vorrichtung nicht in jedem Fall gesichert ist, da die 
Bestimmung der Koordinaten durch Rauschen und Artefakte gestort sein kann. Die 
Position der Vorrichtung wird aus der Differenz von An- und Ausprojektion durch den 

10 Abtastwert mit maximaler Signalamplitude bestimmt. Die Signalqualitat wird jedoch 
durch verschiedene Effekte negativ beeinflusst. Zum einen hangt die Qualitat des 
Signals stark von der Entfernung der Empfangsspule zur Markierung auf der 
interventionellen Vorrichtung ab, da der Impuls uraso schwacher ist, je weiter die 
Empfangsspule vom Signalursprung entfernt ist Einen weit grofieren Einfluss auf die 

1 5 Signalqualitat hat jedoch die Orientierung der Vorrichtung zur Sende- als auch zur 
Empfangsspule. Bei groBen Winkeln zwischen der lokal durch ein Dipolmoment 
angenaherten Schwingspule und den Feldlinien der Sende- bzw. Empfangsspule 
koppeln diese nur noch schwach. 

20 AuBer durch die starke Variation im durch die interventionelle Vorrichtung hervor- 

gerufenen Impuls ist die Lokalisierung wesentlich durch ausgedehnte Aktefakte gestort. 
Haufig wird das Hmtergrundsigna! in der Differenz nicht vollstandig ausgeloscht, was 
darauf zuruckzufuhren ist, dass die Magnetisierung zum Zeitpunkt der Anregung fur die 
jeweiligen Projektionen sich ^nicht ^im gleichen Zustand befindet, sondem einem Ein- 

25 schwingvorgang unterworfen ist Aus diesem Grund sind die Amplituden in der An- 
und Ausprojektion unterschiedlich hoch. Die so hervorgerufenen Artefakte werden im 
folgenden als Einschwingartefakte bezeichnet. 

Weitere Artefakte, die auch als Bildschichtartefekte bezeichnet werden, entstehen, weil 
30 bei jeder neuen Detektion die Magnetisierung in der vorherigen Bildschicht in der Regel 
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nicht voHstandig abgeklungen ist. Diese Restmagnetisierung klingt dann zwischen An- 
nuel Ausprojektion ab und etscheint deshalb in der Dififerenzprojektion als Artetakt in 
der Mitte des Datenvektors. SchlieBUch konnen auch Bewegungen durch Atmung und 
Herzschlag sowie pulsativer Blutfluss negativen Einfluss auf die Quahtat des Signals 
haben. 

Eine gesicherte Aussage fiber die Position der interventioneUen Vorrichtungen ist nicht 
mehr mdglich, wenn sich der durch Rauschen und Artefakte hervorgerufene 
Hintergrund der Amplitude des von der Markierung der interventioneUen Vorrichtung 
ausgehenden Impulses annahert Ausgehend von diesem Stand der Technik stellt sich 
daher die Aufgabe, ein Magnetrescnanzverfehren zur Lokalisierung von 
interventioneUen Vorrichtung zu schaffen, bei dem Rauschen und Artefakte soweit 
unterdruckt werden, dass die Detektierbarkeit des von der Markierung der: 
interventioneUen Vorrichtung ausgehenden Signals stets gewahrleistet ist. 

Die Aufgabe wird erfindungsgemaB gelost durch ein Magnetresonanzverfahren ent- 
sprechend dem Oberbegriff von Anspruch 1, bei dem die Messsignale zur Verbesserung 
der Lokalisierung mit Hilfe eines emdimensionalen Signalverarbeitungsverfahrens 
bearbeitet werden. 

Des weiteren betrifft die Erfindung auch ein Gerat sowie ein Computerprogramm zur 
Durchfuhrung des erfindungsgemaBen Verfahrens. 

Als interventioneUe Vorrichtung werden im Rahmen dieser Erfindung insbesondere 
Katheter, aber auch Biopsienadeln, minimalinvasive chirurgische Instrumente, 
FuhrungsdraTite, Stents usw. verstanden. Die Markierung an der interventioneUen 
Vorrichtung kann insbesondere ein Schwingkreis an der Spitze eines OptiMa-Katheters 
sein; es kann sich jedoch durchaus auch um andere Arten von Anordnungen wie z. B. 
eine Mikrospule handeln, wie sie fur aktive LokaUsierungsverfehren verwendet wird. 
VorteUhaft isthierbei eine ein- bzw. ausschaltbare Markierung, die die separate 
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Aufhahme von Messsignalen im An- bzw. Auszustand erlaubt, im Rahmen dieser 
Erfindung auch als An- bzw. Ausprojektion bezeichnet, so dass die 
Positionsbestimmung der Markierung durch DifFerenzbildung zwischen An- und 
Ausprojektion moglich ist 

5 

Bei dem eindimensionalen Signalverarbeitungsverfahten handelt es sich vorzugsweise 
urn ein iteratives Veriahren, wie es sich fur analytisch nicht direkt losbare Probleme 
anbietet Besonders geeignet ist die sogenannte Maximum-Entropie-Mefhode. 

10 Bei der Maximum-Entropie-Methode (ME-Methode) handelt es sich urn ein iteratives, 
nichtlineares Verfahren zur Signalrestauration. Die ME-Methode lost unterbestimmte 
Probleme, indem sie aus alien mit den Daten vertraglichen Losungen die mit der 
maximalen Entropie wahlt Einen besonderen Vorteil bietet die Moglichkeit, Vorwissen 
fiber den Messprozess durch Einarbeitung zusatzlicher Parameter im Algorithmus zu 

15 berucksichtigen. 

Das der Maximum-Entropie-Methode zugrundeliegende Ausgangsproblem lasst sich 
allgemein wie folgt beschreiben: 

20 Aufgabe ist, eine Verteilungsfiinktion als beste Schatzung jRir eine Verteilung von 
Zustanden zu bestimmen. Oblicherweise gibt es eine unendliche Menge von 
Verteilungen, die mit den Nebenbedingungen vertraglich sind. Das Prinzip der 
maximalen Entropie besagt, dass unter diesen die Verteilung zu wahlen ist, welche die 
maximale Entropie besitzt Diese Auswahl ist die einzig mit den Daten konsistente, 

25 ohne zusatzliche Informationen hinzuzufiigen. 

Ein wahrscheinlichkeitstheoretischer Ansatz zur Begriindung der ME-Methode wird u. 
a. vom G.J. Daniell und S.F. Gull inlEEProc. 127, Pt. E, 170-172 (1980)beschrieben. 
Demnach gilt, wenn das Eingangssignal durch weifies Rauschen uberlagert ist: 

30 
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2 _ — i (Abweichung zwischen MeBsignai und vorhergesagtem Signal) 2 
x " " (Fehler in gemessenem Signal) 2 



Die Wahrscheinlichkeit fur das geschatzte Signal ist dann proportional zu exp(- 
1/2 5c 2 ). Die ME-Mefhode basiert somit auf einer x 2 -Minimierung unter Anpassung der 
geschatzten Signale an die Messdaten. Als besonders geeignet zur Anwendung bei dem 
erfindungsgemaBen Verfahren erwies sich der Algorithmic, der auf die Autoren 
Skilling und Bryan, Mon. Not astr. Soc. 221, 1 1 1-124 (1984) zuruckgeht und sich 
durch eine hohe Konvergenzgeschwindigkeit auszeichnet. 



10 GemaB einer vorteilhaflen Ausgestaltung der Erfindung werden zur Unterdruckung von 
in den Messsignalen auftretenden Artefakten Modellftmktionen gebildet, angepasst und 
bei der Durchfuhrung des iterativen Verfahrens von den Messsignalen subtrahiert Die 
Anpassung der Modellftmktionen an die aufgenommenen Messsignale (die Anprojek- 
tion) erfolgt zweckmaBigerweise, indem die Modellfunktionen mit einem Skalierungs- 

15 parameter angesetzt werden. Die Einbindung in den Maximum-Entropie-Algorithmus 
kann auf zwei verschiedene Weisen erfolgen. Der Skalierungsparameter kann nach 
jedem Iterationsschritt neu angepasst werden oder aber nur einmal vor der ME- 
Interation. In hierzu durchgefuhrten Test wurde im ersten Fall der Parameter in 
Abhangi^ceit des Rauschens mit einer Genauigkeit von 1 bis 4 % bestimmt, wahrend 

20 im zweiten Fall die relative Abweichung ungefahr doppelt so groB war. Auf der anderen 
Seite wird im zweiten Fall ca. 10% weniger Rechenzeit benotigt. 

Einschwingartefekten undBildschichtartefakten unterschieden werden. Da das 
Auftreten von Einschwingartefekten daiauf zuruckzufuhren ist, dass die Magnetisierung 
25 zum Zeitpunkt der Anregung sich fur Messungen, bei denen die Markierung an der 
interventionellen Vorrichtung ein- bzw. ausgeschaltet ist, nicht im gleichen Zustand 
befindet, insbesondere bei Verwendimg des oben beschriebenen OptiMa-Katheters, der 
uber eine optisch schaltbare Markierung verfiigt, wird das Hintergrundsignal durch 
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Bildung der Differenz von Messungen mit aktivierter bzw. desaktivierter Markierung 
nicht vollstandig ausgeloscht. 

Aus diesem Grund bietet sich eine aufgenommene Ausprojektion als ModeUftinktion 
5 zur Unterdriickung der Einschwingartefakte an. Uber den bereits oben erwahnten 
Skalierungsparameter kann die auf diese Weise aufgestellte Modellfunktion an die 
aufgenommenen Messsignale angepasst werden, indem An- und Ausprojektion einander 
gegenubergestellt werden. Wahrend der % 2 -Anpassung wird die Modellfiinktion dann 
vom Messsignal subtrahiert. Das die Position der interventionellen Voirichtung 
10 definierende Signal wird somit relativ zum Hintergrund verstarkt, so dass mit deutlich 
erhohter Sicherheit der Abtastwert mit maximaler Signalamplitude der Position 
zugeordnet werden kann. 

Zur Unterdriickung der ebenfells moglichen Bildschichtartefekte, die darauf zuruck- 
15 zufuhren sind, dass bei den einzelnen Detektionen die Magnetisierung in der vorherigen 
Bildscbicht in der Regel nicht vollstandig abgeklungen ist, miissen hingegen andere 
Modellfimktionen verwendet werden. EBerbei bieten sich Rechteck- oder GauBfunk- 
tionen an, die ebenfalls uber einen Skalierungsparameter angepasst werden konnen. Der 
Grund fiir die Art der verwendeten Modellfunktionen ist im deutlich schmaleren Er- 
20 scheinungsbild der Bildschichtartefakte, die in der GroBenordnung der Bildschichfbreite 
liegen, im Vergleich zu Einschwingartefakten zu sehen. 

Um eine Aussage uber die LeistungsShigkeit der Signalverarbeitung relativ zur Giite 
des Eingqngssignals treffen zu konnen, werden zwei verschiedene Grofien verwendet. 

25 Zum einen liefert das Signal-ZRauschleistungsverhaltnis S/N eine Aussage uber die 
Rauschminimierung nach der Signalverarbeitung, wobei jedoch keine Signalstorungen 
durch Artefakte beriicksichtigt werden, welche unter Umstanden die Bestimmung der 
Position der interventionellen Vorrichtung weit mehr beeintrachtigen als das Rauschen 
alleine. GroBere Aussagekraft hat deshalb das SignaWStorleistungsverhaltnis S/A, 

30 welches zusatzlich zur hochfcequenten Rauschleistung die niederfrequente 



-9- 



PHDE030132 EP-P 



AttefakHeistung mitberiicksichtigt: Hierbei handelt es sich urn den Quotienten aus der 
Nutzsignalleistung und der urn die Leistung des Gleichsignals reduzierten 
Gesamtleistung. Wenn das Rauschen in einem Signal dominant ist, so strebt das S/A- 
gegen das S/N-Verhaltnis. Die Unterdruckung des Rauschens allein hat jedoch nur eine 
5 geringp Verbesserung des S/A-Verhaltnisses zur Folge. Das S/A-Verhaltnis ist deutlich 
besser geeignet als das S/N-Verhaltnis, die Sicherheit der Positionsbestimmung zu 
beurteilen. So hat sich bei den duxchgefiihrten Untersuchungen ergeben, dass eine 
sichere Detektierbarkeit der Position der interventionellen Vorrichtung in aller Regel 
gegeben ist, wenn ein S/A-Verhaltnis >20dB gemessen wird. 

10 

Die Konvergenzgeschwindigkeit des Maximum-Enlropie-Algorithmus ist in erster Linie 
vom Rauschen abhangig. Unabhangig davon kann die Anzahl der Iterationen dvtrch 
geeignete Wahl des benutzerdefinierten Hintergrunds, d. h. des Startwertes der Iteration, 
beeinflusst werden, da je nach Wahl dieses Startwertes die x 2 -Anpassung zu Beginn der 
15 Iteration unterschiedlich gut ist. Eine Erhohung der Konvergenzgeschwindigkeit ist 
besonders wichtig, wenn eine Signalverarbeitung in Echtzeit angestrebt wird. 

Es hat sich herausgestellt, dass bei dem erfindungsgemaBen Verfehren ohne zusatzliche 
Anwendung von Modellfunktionen die Konvergenzgeschwindigjceit maximal ist, wenn 

20 als Startwert der Iteration der Mittelwert der Messsignale gewahlt wird. Gleichzeitig 
ergibt sich fur diese Wahl des benutzerdefinierten Hintergrunds auch das maximale 
S/N-Verhaltnis, wahrend das S/A-Verhaltnis von der Wahl des Startwertes der Iteration 
weitestgehend unabhangig ist Bei geeigneter Wahl des Startwerts konvergiert der ME- 
Algorithmus in weniger als zehn Iterationsschritten. Werden hingegen 

25 Modellfunktionen enteprechend der obigen Ausfiihrungen zur Optimierung der 
Signalverarbeitung und Elimimerung von Artefakten verwendet, hat es sich als 
zweckmaBig herausgestellt, als Startwert der Iteration den Mittelwert der Differenz aus 
Messsignalen und Modellfbnktion zu verwenden. Dieser Mittelwert ist deutlich geringer 
als der Mittelwert des Messsignals, weil die wesentlichen Artefakte bereits durch die 

30 Modellfunktion unterdruckt wurden. 
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Eine weitere Moglichkeit zur Erhohung der Qualitat der Messsignale, die die 
Maximum-Entropie-Methode bietet, besteht darin, Rauschen und Artefakte durch 
Loschen der entsprechenden hoch- bzw. niederfrequenten Eingangssignalanteile zu 
unterdrucken. Da eine sichere Bestimmung der Position der interventionellen 
5 Vomchtung bei Auftreten ausgedehnter Artefakte hoher Amplitude starker 

beeintrachtigt wird als durch das Rauschen allein, ist es besonders wichtig, diese 
Artefakte zu unterdrucken. Sowohl in vitro als auch in vivo wurden vornehmlich 
Artefakte beobachtet, die vier bis fiinftnal breiter waren als der von der Markierung 
ausgehende Impuls. Bei einer Gesamtzahl N von 256 Abtastwerten handelt es sich dabei 
10 typischerweise urn Artefakte, die sich uber mehr als 32 Abtastwerte erstrecken. 

Die Unterdriickung unnotig vieler Signalanteile bewirkt allerdings starke EinbuBen im 
S/N-Verhaltnis, was darauf zuriickzufuhren ist, dass durch Loschen dieser nieder- 
frequenten Signalanteile der Mittelwert stark verringert wird, wahrend das Rauschen im 

15 wesentlichen unbeeinflusst bleibt Entsprechend ist z. B. bei einer Artefakfbreite von 32 
Abtastwerten das S/A-Verhaltnis maximal, wenn 8 niederfrequente Abtastwerte 
eliminiert werden, was dem Quotienten aus der Gesamtzahl der Abtastwerte und der 
Anzahl der Abtastwerte, uber die sich ein Artefokt erstreckt, entspricht Dariiber hinaus 
kann die Loschung zu vieler niederfirequenter Signalanteile, die bei Vorkommen 

20 massiver Artefakte viel Signalleistung enthalten, bewirken, dass die 

Konvergenzkriterien fur den ME-Algorithmus nicht mehr erfullt sind, wenn ein zu 
niedriger Startwert fur die Iteration verwendet wird. 

King VerhesgeriiTi P ; Apt Si p rnwlmialitat <\\vroh Flimirtmmn o r Am PflngoTipmg rmj gntnit V^r-, 

25 besserung des S/N-Verhaltnisses kann durch Loschung von hochflequenten Abtast- 
werten im Spektrum erzielt werden. Die Loschung zu vieler hochfrequenter Abtastwerte 
fuhrt allerdings zu einer starken Verringerung der Nutzsignalleistung, was mit EinbuBen 
im S/A-Verhaltnis verbunden ist. Bei einer Gesamtzahl von N = 256 Abtastwerten 
wurde dabei gefunden, dass nicht mehr als 96 hochfrequente Abtastwerte geloscht 

30 werden sollten, weil in diesem Bereich das Spektrum des Nutzsignals vernachlassigbar 
ist. Ein wesentlicher Einfluss auf die Anzahl der Iterationsschritte durch Unterdriickung 
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hoch- Oder niederfrequenter Signalanteile und damit auf die Rechenzeit konnte nicht 
festgestellt werden. 

In in vivo-Versuchen konnte gezeigt werden, dass durch Eliminierung von Signalteilen 
5 auch bei Vorliegen von stark verrauschten und durch Artefekte gestorten Eingangs- 
signalen eine sichere Positionsbestimmung moglich wird. Es muss jedoch darauf hinge- 
wiesen werden, dass bei der oben beschriebenen erweiterten ME-Methode, in deren 
Rahmen angepasste Modellfunktionen von den Messsignalen subtrahiert werden, die 
Eliminierung von Abtastwerten nicht sinnvoll ist. Dies ist darauf zuruckzufuhren, dass 
10 im Laufe der x 2 -Anpassung die Artefekte entsprechend der Modellfiinktion vom 

Messsignal subtrahiert werden, wobei das Schatzsignal mit diesem Differenzsignal in 
Ubereinstimmung gebracht werden muss. Eine zusatzliche Loschung niederfrequenter 
Signalanteile wurde daher zu einer Verfalschung fuhren, die eine Anpassung nicht mehr 
erlaubt 

15 

Neben den oben beschriebenen iterativen Verfehren, insbesondere der Maximum- 
Entropie-Methode, kdnnen auch andere eindimensionale Signalverarbeitungsverfahren 
wie z. B. Filter zum Einsatz kommen. Grundsatzlich geeignet sind sowohl Filter mit 
endlicher als auch solche mit unendlicher hnpulsantwort, die man auch mit den 
20 englischen Begriffen FIR (finite impulse response) bzw. IIR (infinite impulse response) 
bezeichnet. Derartige Filter sind dem Fachmann grundsatzlich bekannt. Zwei typische 
Filter, die sich als geeignet zur Losung der erfindungsgemaQen Aufgabe herausgestellt 
haben, sind das Wiener-Filter und das Bandpassfilter. 

25 Das Wiener-Filter kann in Fourierform dargestellt werden als: 
W 1. ggtf 
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Dabei sind H die tTbertragungsfunktion des MeBsystems und O ff und 3>nn die Leistungs- 
dichtespektren des gesuchten Signals fk bzw. des Rauschens nk. 

Das Wiener-Filter ist besonders geeignet, das S/N-Verhaltnis zu vefbessern, d. h. das 
5 Rauschen effektiv zu unterdriicken. Artefakte hingegen werden weniger gut unterdruckt 
als bei Anwendung der Maximum-Entropie-Methode . 

Bin weiteres geeignetes Filter ist das Bandpassfilter, das sich als wirkungsvoll zur 
Unterdruckung von Rauschen und Artefakten erwiesen hat. Die Sicherheit bei der 
10 Lokalisierung einer interventionellen Vorrichtung konnte mit Hilfe eines 

Bandpassfilters deutlich gesteigert werden. Lediglich zur Unterdruckung schmaler 
Artefakte, wie z. B. Bildschichtartefakte, ist der Bandpassfilter weniger gut geeignet. 

Die Wahl des am besten geeigneten Signalverarbeitungsverfahrens ist abhangig von der 
15 exakten Problemstellung. Auf der einen Seite zeigt die Maximum-Entropie-Methode 
die besten Ergebnisse hinsichtlich der Artefakt- und Rauschunterdriickung, besonders 
bei Implementierung der zusatzlichen, oben ausfuhrlich dargelegten Merkmale. 
Andererseits erfordert die ME-Methode als iteratives Verfahren deutlich mehr 
Rechenzeit als die Anwendung eines Filters. Wahrend diese fur die Filter im Bereich 
20 von 1 bis 2 ms liegt, kann die Rechenzeit je nach Gesamtzahl der Abtastwerte fur die 
ME-Methode > 100 ms sein. Bei sehr hohen Anforderungen hinsichtlich der Kurze der 
Rechenzeit fur eine Echtzeitvisualisierung kann daher ein Filter der ME-Methode 
vorzuziehen sein. 

~25 Eine weitere Verbesserung bei der Lokalisierung einer interventionellen Vorrichtung 
kann bei Vorliegen mehrerer, der Lokalisierung dienender Messsignale dadurch erzielt 
werden, dass nach Bearbeitung der Messsignale mit Hilfe des eindimensionalen 
Signalverarbeitungsverfahrens eine Uberpriifung der uber die bearbeiteten Messsignale 
ermittelten Positionen der interventionellen Vorrichtung auf Koinzidenz durchgefiihrt 
30 wird. Eine solche Uberprufung bietet sich insbesondere bei Verwendung des oben 
beschriebenen OptiMa-Katheters an, bei dessen Verwendung sich auf dem Korper des 
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Patienten mehrere Empfangsspulen befinden, die die Messsignale parallel aufnehmen. 
Diese Messsignale unterscheiden sich bei der Lokalisierung zwar hinsichtiich der 
Amplitude, sollten jedoch die gleiche Position im Raum fur die mterventionelle 
Vorrichtung ergeben. 

5 Bei einer Uberpriifung der bearbeiteten Messsignale auf Koinzidenz wird nun im 
Anschluss an die Bearbeitung der Messsignale gepruft, ob die fiber die einzelnen 
Empfengsspulen ermittelten Positionen ubereinstimmen. Eine solche vollstandige oder 
teilweise Ubereinstimmung erhoht die Wahrscheinhcbkeit, dass die ennittelte Position 

10 korrektist, zusatzlich. 

Vorzugsweise erfolgt die Verarbeitung der verschiedenen, der Lokalisierung der 
interventionellen Vorrichtung dienenden Messsignale gemeinsam in dem 
eindimensionalen Signalverarbeitungsverfahren, so dass auch die Einflusse auf die 

15 Positionsbestimmung fur die einzelnen Messsignale gleich sind. Dies ist sowohl durch 
Anwendung ernes iterativen Verfabrens wie der Maximum-Entropie-Methode als auch 
durch Verwendung eines Filters mogUch. AnschlieBend kann die tTberpriifung der 
ermittelten Positionen fur die interventionelle Vorrichtung auf Koinzidenz durchgefuhrt 
werden. Es kann auch direkt durch das emdimensionale Signalverarbeitungsverfahren 

20 die Korrelation der Messsignale berechnet werden, urn auf diese Weise ein MaB fur die 
Obereinstimmung der Signalspektren zu erhalten. 

Die Erfindung wird anhand der beigefugten Figuren weiter verdeutiicht. Es zeigen: 

25 Fig. 1. Auftragung der Signalamplituden gegen die Abtastwerte zur Veranschaulichung 
der Signahestauration mit Hilfe der erweiterten ME-Methode bei starker 
Stoning der Eingangssignale durch Emschwingartefakte; 



Fig. 2 : 

30 



Auftragung der Signalamplituden gegen die Abtastwerte zur Veranschaulichung 
der Signalrestauration mit Hilfe der erweiterten ME-Methode bei starker Stoning 
der Eingangssignale durch Bildschichtartefakte. 
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In Figur 1 (a) ist ein in vitro-Eingangssignal mit einer Gesamtzahl von N = 256 
Abtastwerten gezeigt, in dem die Katheterposition durch einen Pfeil markiert ist. Die 
Signalamplituden auf der Ordinate sind fur die einzelnen Abtastwerte auf der Abszisse 
5 graphisch dargestellt Die Messungen warden mit Hilfe eines 1,5 Tesla MR-Tomo- 
graphen (GyroScan ACS-NT, Philips Medical Systems) unter Anwendung einer 
„spoiled" Gradient-Echo-Sequenz (FOV = 256 mm) durchgefiihrt, wobei der Katheter, 
bei dem es sich urn einen OptiMa-Katheter handelte, in einem Schlauchphantom 
platziert wurde. Das Eingangssignal ist stark durch Einschwingartefakte gestort, die 

1 0 durch Bildung und Anpassung einer Modellfunktion eliminiert werden, welche im 
Rahmen des ME-Verfahrens von den Messsignalen subtrahiert wird. Als 
ModeUfunktion wird die in (b) dargestellte Auspiojektion verwendet, bei der es sich urn 
die aufgenommenen Signale bei desaktivierter Markierung am Katheter handelt In (c) 
ist das Ergebnis nach erfolgreicher Signalrestauration dargestellt, woraus die eindeutige 

15 Bestimmbarkeit der Katheterposition deutlich ersichtlich ist Die Signalverarbeitung ist 
mit einem erheblichen Anstieg des S/N- wie des S/A-Verhaltnisses verbunden. In 
analoger Weise wild in den Darstellungen (d) - (f) die Signalrestauration eines stark 
durch Einschwingartefakte gestorten in vivo-Eingangssignals gezeigt, wobei hier die 
Gesamtzahl der Abtastwerte bei N = 128 lag. Fig. 1 (d) zeigt hier das Eingangssignal, 

20 (e) die entsprechende Ausprojektion und (f) das Ergebnis nach erfolgter 

Signalrestauration. Fur die in vivo-Messungen wurde das gjeiche Verfahren wie fiir die 
in vitro-Messungen angewandt, wobei hier jedoch ein entspiechender Katheter in die 
Aorta eines Schweines eingefiihrt und eine refokussi erte Gradient-Echo-Sequenz (FOV 
= 300 mm) verwendet wurde. 

25 

Figur 2 (a) zeigt ein Kathetersignal mit schmalen Bildschichtartefakten, wobei 
wiederum die Signalamplituden fiir die einzelnen Abtastwerte graphisch dargestellt sind 
und die Position des Katheters mit einem Pfeil gekennzeichnet ist Die im Rahmen der 
erweiterten ME-Methode verwendeten Modellfunktionea, die bei dem iterativen 
30 Verfahren wiederum von Messsignalen subtrahiert werden, sind in (b) dargestellt In (c) 
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ist zu sehen, dass nach Signalrestauration die Position des Katheters eindeutig bestimmt 
werden kann, obgleich die in (a) auftretenden Artefakte sehr schmal sind und die wahre 
Katheterposition in der Amplitude ubertreffen. 
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PATENTANSPROCHE 



1 . Magnetresonanzverfahren zur Lokalisierung von interventionellen Vorrichtungen, 
insbesondere in vivo, bei dem die interventionelle Vomchtung eine Markierung tragt, 
die bei Magnetresonanzaufoahmen die Messsignale beeinflusst oder eigene Messsignale 
5 erzeugt, 

dadurch geken ngeiclitiet, 

dass die Messsignale mit Hfflfe eines eindimensionalen Signalverarbeitungsverfahrens 
bearbeitet werden. 

10 2. Verfahren nach Anspruch 1 , 
dadurch geket mrairhTiet, 

dass das eindimensionale Signalverarbeitungsverfehren ein iteiatives Verfahren ist 

3. Verfahren nach Anspruch 2, 
15 dadurch gekennzeichnet 

dass das iterative Verfahren auf der Maximum-Entropie-Methode basiert. 

A V^rfbVirP!n n»rVh Anwrnoli 0 rw^r^ - : — : : : 

, _ , f — — ■ — 2" " n '. 9 

dadurch geketi ti^eichnet, 
20 dass fiir in den Messsignalen auftretende Artefakte Modellfunktionen gebildet, 

angepasst und bei der Durchfiihrung des iterativen Verfahrens von den Messsignalen 
subtrahiert werden. 
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5. Verfahren nach Anspruch 4, 
A*An*fh preVftrniseiclmet. 

dass die Modellfunktionen iiber einen Skalierungsparameter an die aufgenommenen 
Messsignale angepasst werden. 

5 

6. Verfahren nach Anspruch 5, 
dadjggh e& trnnveichneL 

dass die Anpassung der Modellfiinktionen an die aufgenommenen Messsignale nach 
jedem Iterationsschritt im Rahmen des iterativen Verfahrens neu erfolgt 

10 

7. Verfahren nach Anspruch 5, 
ria/inTr* pwVftnr>7:eichnet. 

dass die Anpassung der Modellfunktionen an die aufgenommenen Messsignale einmalig 
vor der Durchfuhrung des iterativen Verfahrens erfolgt. 

15 

8. Verfehren nach einem der Anspruche 4 bis 7, 

dadUTCh g j^remn^chneiL 

dass als Modellfunktion die bei inaktiver Markierung an der interventionellen 
Vorrichtung aufgenommenen Messsignale verwendet werden. 

20 

9. Verfehren nach einem der Anspruche 4 bis 8, 

dass als Modellfunktionen Rechteck- oder GauBfunktionen verwendet werden. 

25 10. Verfehren nach einem der Anspruche 4 bis 9, 
dadurgh pekerm geichnet. 

dass als Startwert der Iteration der Mittelwert der Differenz aus Messsignal und 
Modellfunktion gewahlt wird. 
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1 1 . Verfahren nach einem der Anspriiche 2 bis 9, 
dadurch flekerin geichnet 

dass als Startwert der Iteration der Mittelwert des Messsignals gewahlt wird. 
5 12. Verfahren nach einem der Anspriiche 1 bis 1 1 , 

riariiirch p;ekenti7 :eichnet. 

dass zur Unterdruckung von Rauschen und/oder Artefakten in den aufgenommenen 
Messsignalen hoch- und/oder nieder&equente Signalanteile elirniniert werden. 

10 13- Verfahren nach Anspruch 1, 
dadurch gekanti ^eichnet 

dass als eindimensionales Signalverarbeitungsverfahren ein Filter mit endlicher oder 
unendlicher Impulsantwort verwendet wird 

15 14. Verfahren nach Anspruch 13, 
dadurch gekennzeichnet 

dass das Filter ein Wiener-Filter oder ein Bandpassfilter ist 

15. Verfahren nach einem der Anspriiche 1 bis 14, 
20 dadurch gekennzeichnet 

dass bei der Auswertung mehrerer, der Lokalisierung der interventionellen Vorrichtung 
dienender Messsignale nach Bearbeitung der Messsignale mit Hilfe des 
eindimensionalen Signalverarbeitungsverfahrens eine tJberpriifung der uber die 
bearbeiteten Messsignale ermittelten Positionen der interventionellen Vorrichtung auf 

ry tr vr - " • i . i i r~—* . • ■+ : 

16. Verfahren nach einem der Anspriiche 1 bis 15, 
dadurch geke tinrairrtiiiftt, 

dass mehrere, der Lokalisierung der interventionellen Vorrichtung dienende 
30 Messsignale gemeinsam im eindimensionalen Signalverarbeitungsverfahren bearbeitet 
werden. 
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17. Verfahxen nach einem der Anspruche 1 bis 16, 
dadurch gp Venngeichnet. 

dass die Messsignale von mehreren Empfangsspulen parallel aufgenommen werden. 

18. Verfahren nach einem der Anspruche 1 bis 17, 
dadurch p^Vgmirafohnet. 

dass das eindimensionale Signalverarbeitungsverfahren die Korrelation einer oder 
mehrerer Messsignale berechnet 



19. Gerat zur Lokalisierung von interventioneUen Vorrichtungen mit Hilfe von 
Magnetresonanzaufnahmen, bei dem die interventionelle Vorrichtung eine Markierung 
trSgt, die bei den Magnetresonanzaumahmen die Messsignale beeinflusst oder eigene 
Messsignale erzeugt, 

15 dadurch ge VftnnTOiclmet 

dass das Gerat fiber eine Programmsteuerung zur Durchtuhrung eines Verfehrens nach 
einem der Anspruche 1 bis 18 verfugt. 

20. Computerprogramm zur Verarbeitung von Messsignalen bei der Lokalisierung von 
20 interventioneUen Vorrichtungen mit Hilfe von Magnetresonanzaumahmen, bei dem die 

interventionelle Vorrichtung eine Markierung trSgt, die die Messsignale bei den 
Magnetresonanzaumahmen beeinflusst oder eigene Messsignale erzeugt, 
fa^yrrfo prftVflttnr eichnet 

dass mit dem Computerprogramm ein Verfahren nach einem der Anspruche 1 bis 18 
25 durchfuhrbar ist. 
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ZUS AMMENF AS SUNG 
Magnetresonanz-I^kaUsierimgsverfahren 

Die Erfindung betrifft ein Magnetresonanzverfahren zur Lokalisierung von interven- 
tionellen Vorrichtungen, insbesondere in vivo, bei dem die interventionelle Vorrichtung 
5 eine Markierung tragt, die bei Magnetresonanzaufoahmen die Messsignale beeinflusst 
oder eigene Messsignale erzeugt, wobei zur Unterdriickung von Rauschen und Arte- 
fakten eine Verarbeitung der Messsignale mit Hilfe eines eindimensionalen Signalverar- 
beitungsverfahrens erfolgt. Hierbei kann es sich insbesondere urn die Maximum- 
Entropie-Methode bandeln, die durch Verwendung von Modellfunktionen noch erwei- 
10 terbar ist Diese Modellftmktionen werden im Rahmen des iterativen Verfahrens von 
den Messsignalen subtrahiert, um auf diese Weise die Eliminierung von Artefekten 
zusatzlich zu verbessern. Alternativ zur Anwendung der Maximum-Entropie-Mefhode 
ist auch die Verwendung von Filtern, insbesondere Wiener- oder Bandpassfiltern 
moglich. 

15 

Fig. 1 
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